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Abstract 

We consider the dynamics of cargo driven by a collection of interacting molecular motors 
in the context of an asymmetric simple exclusion processes (ASEP). The model is formulated 
to account for i) excluded volume interactions, ii) the observed asymmetry of the stochastic 
movement of individual motors and iii) interactions between motors and cargo. Items (i) and 
(ii) form the basis of ASEP models and have already been considered in the literature to study 
the behavior of motor density profile [1]. Item (iii) is new. It is introduced here as an attempt 
to describe explicitly the dependence of cargo movement on the dynamics of motors. The 
steady-state solutions of the model indicate that the system undergoes a phase transition of 
condensation type as the motor density varies. We study the consequences of this transition 
to the properties of cargo velocity. 
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1 Introduction 



Asymmetric simple exclusion processes (ASEP) are specially convenient for describing general prop- 
erties of dynamical systems consisting on a collection of many-interacting particles in situations for 
which the physicochemical characteristics of the components and thus the nature of interactions, 
need not to be described in details [1], [2]. Because of this, ASEP models have been used to study 
the collective movement of molecular motors that happen at the microtubules within cellular en- 
vironment [3], [1], [5], [6]. These are models that can be defined in one-(spacial)-dimension and 
incorporate the asymmetry of the motion of individual motors. 

Since the initial proposal pointing out to ASEP as a possibility to describe the collective dy- 
namics of molecular motors [3], the general interests are mainly focused on the properties of the 
system at different boundary conditions that allow to make predictions on the stationary currents 
or on the average motor velocities as function of external loading forces [7]. Also, the effects of 
motor coordination onto the process of pulling on fluid membranes ^ have been studied in the 
literature in the context of discrete ASEP models with disorder [9], ^U\. It shall be interesting then 
to investigate whether this type of description can be extended to describe directly the movement 
of cargo driven by motors to understand some of its characteristics observed in experiments. 

Cargo transport by motors happens at the cellular environment where simple diffusion of vesicles 
or nutrients is severely limited by the presence of innumerous structures inside the cytoplasm. 
Besides, it is known that virus particles can take advantage of the existing transport mechanisms 
using molecular motors as carriers to reach the interior of the cell [TT], [12]. Therefore, by studying 
the properties of a model that leads to quantitative predictions on the movement of both motor 
and cargo might be helpful to understand mechanisms to prevent cell infection and/or to design 
more efficient drug carriers |13j . 

We have already worked on this problem to investigate the movement of cargo in connection to 
the short-time behavior of the motor density profile defined in the context of the continuum limit 
of an ASEP with periodic boundary conditions [H]. Here, we present an alternative to describe the 
long-time regime (steady-state) of the movement of cargo using a discrete version of the model. The 
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elementary dynamical processes, that is, the processes at the level of individual particles moving 
on a defined one-dimensional lattice with periodic boundary conditions are such to account for 
two kinds of particles - the motors and the cargoes. Other processes that take place on a lattice 
and involve two kinds of particles with exclusion have already been explored in the literature to 
describe diverse phenomena. The presence of the seminal second class particle, for example, has 
been considered to study the microscopic properties of shock fronts exhibited by the system due to 
density inhomogeneities [15]. This kind of model can also be used to study the properties of the 
density profile in the presence of defects [T6] . 

The important point we want to notice is that, up to the present, all kinds of particles in 
the ASEP models presented in the literature are, in all cases, provided with their own (intrinsic) 
dynamics. As a rule, an encounter of a particle with another is just supposed to change their 
original dynamics, or even to impede their movement (excluded- volume interactions), by inducing 
modifications on the rates of the original stochastic movement. This is the case, for example of the 
model cited above containing second class particles. The elementary processes in this case are 



10 ^ 01 

12 21 first and second class (1) 
20 02 

that is, either particles of first (1) and second (2) classes can move by their own if the target 
neighboring site is unoccupied (0). They can also move if interacting with other particles by 
interchanging occupancy sites as in 12 ^ 21. This is also the case of another ASEP-like model for 
which the particles with different dynamics simulate the presence of cars (1) and trucks (2) on a 
same traffic road [17J, [18j. For this model, the direction of the intrinsic movement of particle (2) 
is opposed to that in the previous case: 

10 ^ 01 

12 21 cars and trucks (2) 
02 20 

As we want to describe the movement of cargo driven by motors using an ASEP-like model, we 
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need first to think on possible ways cargo may use the motors to move and up to what extent the 
presence of cargo can affect the intrinsic dynamics of motors. The first idea that occurred to us 
is that the elementary processes in this case should be such to attribute exclusively to one kind of 
particles - the motors - the ability to move independently, i.e. motors should be able to hop from 
site to site the only restriction being excluded volume interactions with other motors. The other 
kind of particles - the cargoes - should move only if assisted by motors. Therefore, the dynamics 
must, in addition, incorporate explicitly a form for the motor /cargo interactions at this particle 
level. We want to examine how such interactions affect the long-time average properties of both 
cargoes and motors. 

Based on the above considerations we present in Section 2 what we conceived as a minimum 
model that is able to account explicitly for the dynamics of both cargo and motors expected as a 
result from their mutual interactions. This is presented as an ASEP model that incorporates a few 
common characteristics of this biological system, but avoids details of motor/cargo interactions. It 
turns out that this model is exactly solvable, that is, the steady-state configurations can be deter- 
mined. Using the matrix approach proposed and developed by Derrida, Lebowitz, Evans, among 
others [2], [E], these states are conveniently represented as products of certain non-commuting 
matrices which can be used to calculate the properties of interest as the average cargo velocity 
discussed in Sec. 3. The analytical results obtained in this way are used in Section 4 to discuss the 
phenomenological consequences of the model. 



2 An ASEP model for motors and cargo 

There are a few proposals in the literature to characterize the origins and the role of the components 
involved in the cargo/motor interactions at a microscopic level [20]. Data suggest that such 
interactions are mediated by certain proteins - dynactin is the most studied - but apparently there 
is consensus about their short-range nature as a general characteristic. In terms of the scales 
involved, this is equivalent to say that such cargo/motor interactions happen by "direct contact" 
among individual components. Moreover, although the experiments indicate that the number 
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of motors attached simultaneously to cargo may be large [21], there is limited information on 
the typical times each motor, or group of motors, remains attached to cargo in the course of its 
movement. Therefore, in building up our model, we avoid details of the processes associated to 
such microscopic interactions and simply account for these as stochastic processes. We suppose 
that the time scales associated to the motor/cargo interactions are very short if compared with 
observation time so the events that happen within such intervals need not to be described as we 
treat the problem at a larger time scale. 

The ASEP model considered here consists on a collection of M motors, and K cargoes, referred 
in the following as particles of type 1 and particles of type 2 respectively, distributed among the 
sites of a one dimensional cyclic lattice. Each particle occupies a single site and — M — > 
sites remain empty (Fig.l). The total number of particles of each species is conserved. We consider 
the case A' = 1 so that at the steady state, all configurations of the system are equally likely (p). 
Each site is identified by its position on the lattice j = 1,2,...N and occupation at each site is 
specified by a corresponding site variable aj that assumes integer values 0, 1 or 2, if the site is 
empty (0), occupied by a motor (1), or occupied by cargo (2). A state C of the system is specified 
by the set {(Jicr2... cn}- As usual, in this kind of description the stochastic dynamics is defined 
through a Poissonian process taking place at the lattice such that at each time interval dt, a pair of 
consecutive sites i and i + 1 is selected at random and the system is updated depending on whether 
it is possible to exchange particles between these two sites. We choose the following possibilities 

10 01 with rate k, probability kdt 

12 — s> 21 with rate w, probability wdt (3) 
21 — i> 12 with rate p, probability pdt 

where the pair of sites + is being represented by the values assumed by the variables (cXi, cTj+i). 
According to these rules a cargo is allowed to move only if "assisted" by a motor at a neighbor 
site. We see this as a possibility to describe the fact that the movement of cargo is conditioned to 
that of motors, as observed in real systems. At this level of description one must consider that if 



^In the presence of more than one particle of type 2 obeying the algebra (O, the system loses ergodicity. In these 
more general cases, an invariant measure should be assigned to each subspace of configurations that preserve the 
number of empty sites between each pair of these particles. 



motors affect the movement of cargo then cargo must have influence on the movement of motors. 
In the present model this is incorporated into the dynamics ([3]) both explicitly, by the processes 
that involve cargo/motor exchanging positions in both directions and implicitly, by modifying the 
motors hopping rates that assume distinct values depending whether a jump occurs towards an 
empty site or by interchanging position with particle 2. 

Here, we are concerned with the kinematics of the cargo at long-time regimes. For this, we use 
the matrix-approach [I], [2] to represent any configurations of the system of sites by a product 
of matrices, so that the probability of occurrence of a particular configuration C is given by 

1 ^ 

Pnm{C) = Tr T\{6,^^,D + d^^^^A + 6,^^,E) (4) 

where the Kronecker delta symbols select the correct occupancy and 

AT 

is the normalization. The sum extends over all allowed configurations for which Sa,,i = M and 
X^r^o-i.s = K = 1. In the product above, site i is represented by matrix D if it is occupied by a 
motor (particle 1), by matrix A if it is occupied by the cargo (particle 2) or by matrix E if it is 
empty (not-occupied). 

In the stationary state, the probabilities PN,MiC) for all configurations C satisfy the condition 

Y,Pn,m{C')W{C' ^ C) - Pn,m{C)W{C ^C') = (6) 
c 

where W{C' C) is the rate at which the exchange of particles occur between neighboring sites 
so that all nonzero terms in the above sum are those for which configurations C and C differ from 
each other at most by the occupancy of a pair of sites. 

As an example, we consider for A^ = 4,i^'=l,M = 2 the following configuration C = 1201. 
There is just one way to leave this configuration that is through the process: 12 ^ 21 with 
W{C —>■ C) = w. On the other hand, there are two ways to reach this configuration, either 
by exchanging particle positions (i) by the process 10 01 in configuration C = 1210 with 
WiC ^C)=k or (ii) by the process 21 ^ 12 in configuration C = 2101, with W{C' ^ C) = p. 
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So, in this case, equation ([6]) reads: 

wTr{DAED) = kTr{DADE) + pTr{ADED) 



(7) 



As a general rule, the main difficulty in using this method to determine the probabilities Pn,m{C) 
for each configuration C is to find the algebra, if any, that must be satisfied by the corresponding 
matrices of a given ASEP model in order to satisfy condition ([6]) for a given dynamics as in ([3]). In 
the present case, we conjecture that if D, A and E are such that 

DA - xAD = E-D 
DE = E 

(8) 

EA = E 
EE = E 



for 



X = (9) 

w 



then, Eq. is satisfied. This can be tested using explicit examples, as the one considered in 
(jTj). Using ([HD to evaluate the traces, one can easily check that the identity holds trivially. In the 
following, we study properties of this model that are of interest for examining the consequences of 
the model regarding the characteristics of cargo movement. 



3 Average cargo velocity 



For K = 1 i.e. just one cargo, in the presence of M motors distributed along a cyclic lattice of 
sites, we consider A^ > M + 1, to ensure that at least one site in the system remains empty. In this 
case, the average cargo velocity < f > at steady state is expressed as 

^ / N-2 N-2 \ 

< ^ >= ^ — p 5Z n ("^-^-i^ + ^-..0^)^/^ - ^ 5Z n ("^--1^ + ^-^fiE)DA (10) 

where the sums in the numerator extends over all configurations of M — 1 motors distributed 
among N — 2 sites. The first sum in the RHS accounts for all configurations in which the site at the 



immediate right of cargo is occupied by a motor. The second sum accounts for all configurations 
in which there is a motor at immediate left of cargo at a lattice position. Notice that due to the 
invariance of the trace the normalization factor Zm,m can be written as 

N-l 

= 5^^^ n ("^--1^ + ^-.-0^)^ (11) 

where the sum extends over all configurations of M motors distributed among — 1 sites. 

In order to make reference to the above traces over products of matrices, it is convenient 
to introduce the functions Wa^_^,aj aj+i to indicate the configurations having the n-uple (j — 
+ 1, ...) fixed, the corresponding sites occupied by particles or holes assigned by the variables 
cTj-i, CTj, cTj+i, .... Using these definitions, we write 



— \ pJ^W^i - wJ^Wu] (12) 

V,M \ 7^ s\ I 



<v>=- 
where 

5^W^2i + 5^W^2o (13) 

N-2 

for Wi2 = n (^o-i,i-D + 6^^fiE)DA and analogous definitions for W21, W20 and PFo2- Alternatively, 

i=2 

due to the cyclic property of the trace, Z]\r^M can also be calculated from 

5^W^12 + 5ZVr02. (14) 

These two expressions ( !T3|) and (IT^ are equivalent and both will be used below, at convenience. 
Notice that we can rewrite the sum over configurations in W20 as 

j2w2o = J2 ^120 + J2 ^020 (15) 

Wi} {a-i} {cTi} 

where the first (second) sum on the RHS extends over all configurations of M — 1 motors with 
the triplet 120 fixed (M motors with the triplet 020 fixed) distributed among — 3 lattice sites. 
Making use of the decompositions in f|T3|) and (|T5|) . the average cargo velocity f|T2|) is expressed as 



<V>=p--^\p\J2Wl20 + Y.^02o] + w[j2Wl2]\ (16) 




We compute Zjy^M as it is expressed in Eq. (fH|) . 

A convenient way to perform the calculations indicated above is to replace the sums over site 
variables {o"i} by sums over blocks defined by the integers {rrii} and {qi} for i = l,2...k. (see for 
example, ref.[T^). In this representation, 

J2Wi2= Yl tr{E'''D"'K..E''^-W""=-'E''^D'^>'A) (17) 

with TUk > 1; 

^1^02= tr{E'"D"'\..E'''=-W"'^-'E'''A) (18) 

with gfc > 1; 

^1^120= tr{E'''D'^\..E''^-W"'^-^E'i^D"'^A) (19) 
with gi > 1 and nik > 1; 

YWo2o= Y tr{E'^W"'\..E'^^-W"'^-^E''^A) (20) 

Wi} {mi};{qi} 

with gi > 1 and Q'a,. > 1; 

From the algebra in ([H]), it follows that D"^AE = x^AE. This identity is needed in the evaluation 
of the above traces for general configurations of the variables The results are 

Wi2 = tr{E'iW^\..E'i^-W^^-^E'i^D'^^A) = x'^Hr{E) 
Wq2= tr{EiW^\..Ei>'-W^''-^E'i^A) = tr{E) 
Wi2o = tr{E''W^\..E'i*'-W^^-^E'i^D"'^A) = x"'Hr{E) 
Wo2o= tr{E''W^\..E'i^~^D"'^-^E'i^A) = tr{E) 
Notice that W12 and W120 are functions of the size of the block, i.e. of the number of particles 
of type 1 (motors) that precede particle 2 (cargo). Because of this, the evaluation of the respective 
sums over {qi} and {m,} in configurations of the type Wu (or 1^120), excluding m^, is equivalent 
to account for the number of ways for distributing M — rrik motors into — — 2 (or into 
N — iTLk — 3) sites. The factor 2 in the first case comes from the exclusion of two sites from the 
total: one occupied by the cargo and another that must remain empty to define the limits of the 
cluster of motors behind the cargo. Then, 

(22) 

{a,} m=l \ / 
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and 

M-m 



(I 

{aj m=l ^ 



x™tr(E) (23) 



Wo2 and W020 correspond to configurations that do not present motors behind the cargo. In the 
sum over configurations of the type W02 one must account for the number of ways to distribute M 
motors into N — 2 sites (from the total of sites, there must be excluded 2, one to fix the cargo 
and the other to fix an empty site). Then, 

E^02=f^~^Vr(i?) (24) 



and 



Y.^02,= ,Atr{E) (25) 



We now proceed by computing the sum over integer m in (l22l) and (l23l) . 

3.1 Aproximate expression for the average velocity of the cargo in the 
Umit of very large systems 

Our intention is to obtain an expression for the average velocity < f > of cargo in the limit for 

which both the number of sites and the number of motors (which are conserved by the dynamics) 

are taken very large, that is 00 and M — * 00. These limits are supposed to be taken in such 

a way to ensure that the ratio between these two quantities 

M 

lim — = p (26) 

M^oo 

converges to a defined density of motors p, such that < p < 1. These same limits have already 
been considered in Ref. [T7j to calculate the average velocity of trucks in a related traffic problem. 
Here, we proceed along the same lines sketched by these authors. 



First, we use Stirling formula A^! ~ \/2'kNN^ exp(— A^) for very large A^, to approximate the 
combinatorial coefficients in (I23D and (ESI). It results 



C„. . lu„ i"" - ^\ ~ lim i^-P?-M-N(p^np+(l-p)Ml-m 



N-^QO 

M^oo 
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C„.„ . lim ~ li,„ (l-p)^exp[-iV(,ln,+ (l-p)ln(l-p))| 



and 



In order to evaluate the sums 5*1 and 5*2 in (1221) and (l23l) . we follow the procedure used in Ref . [22] . 
There, sums involving factorials of this kind are approximated by integrals and the asymptotic 
regimes are obtained using Laplace's method [23j. Considering then the limit of very large systems 
and defining z = m/N, we calculate 

lim ^1 ^ lim (-^] 4 (29) 

where we have defined the function f{z) of a single variable z as 



f{z) = (l-2)ln(l-;z)-(p-2)ln(p-;z)-(l-p)ln(l-p) + 2lnx (30) 

M p 

and used the fact that in the specified limit, the sum in m converges to the integral as ^ ^ J dz. 

m 



Now observe that f{z) has a maximum at 

1 — xp 



1 — X 



(31) 



so in order to apply Laplace's method in the present case, one must distinguish between two 
possibihties, namely 

• if xp > 1 (a; > 1) then, < Zma.x < P i-e. ^max belongs to the integration interval. In this 
case, Laplace's method gives 

lim )exp[-Ar(-lnx + (l-p)ln(x-l))] (32) 



N^oo 
M-*oo 



or 



if xp < 1 either for x > 1 or x < 1, then z^g^^ < 0. Therefore, in this case ^max does not 
belong to the integration interval. Since f{z) is a monotone decreasing function of z, the 
integral is dominated by the value of the integrand at z = 0, and application of Laplace's 
method results 

lim S, }\~ f — ^ exp[-iV(p In p + (1 - p) ln(l - p))] (33) 

^2nNp{l-p) In(px) 
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Analogously, the asymptotic behavior of the sum 5*2 (123|1 must be analyzed according to the 
range of xp: 



if xp > 1 {x > 1), then 



lim 5^2 ~ 

iV— ►cxa 
M^oo 



X — 1 



X 



exp[-A^(-lna; + (1 - p) \n{x - 1))] 



(34) 



or 



if xp < 1, (either x > 1 or x < 1), then 



lim 5^2 ~ 

N^oo 



27riVp(l -p) In(px) 



exp[-Ar(plnp+ (1 -p)ln(l -p))] (35) 



From these results, one concludes that if xp > 1 then both 5*1 and 5*2 are the dominant factors 
in the expression for < v > both in the numerator and in the denominator. In this regime. 



lim < V >~ lim 

M^oo M—*oo '- 



kw 
k + p 



for xp > 1 



(36) 



On the other hand, if xp < 1, all factors in the expression (ITB]) for < v > are of the same order of 
magnitude and then, 

< V >~ lim <p— -—— 

M^oo 



lim <v>r^ lim ^ p - [p {S2 + C020) + wSi] ^ ~ pp + 

iV^oo Af^oo I [01 + U02) 



W 



In(px) — 1 



for xp < 1 



(37) 



The consequences of these results to the phenomenology of cargo movement will be analyzed in the 
next section. 



4 Discussion and concluding remarks 

The aim of the present work is to study the consequences of introducing cargo into certain lattice 
models where there is also present a set of biased molecular motors interacting through excluded- 
volume interactions. We consider an ASEP-like model specially designed to account for both 
kinds of particles. Therefore, the model includes an assumption about the stochastic nature of the 
movement of cargo and its dependence on the dynamics of motors. To our knowledge, this is the 
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first attempt to include in the same framework, at the particle level, the effects on the movement 
of motors due to the presence of cargo and vice-versa. 

We look for the probabilities associated to the configurations of the system at the stationary state 
which are represented by products of certain noncomuting matrices pQ. Using this representation, 
we were able to make quantitative predictions on the average properties that characterize the 
movement of cargo. We focus on the computation of the average velocity of cargo < i; > whose 
behavior predicted by the model suggests that the system displays a phase transition under variation 
of the parameters. Fig. 2 shows < v > plotted according to the results in Eqs. (!36l) and (!37|) as a 
function of p, at different values of p for fixed w and k. Observe that for sufficiently high values 
of p the function < v > displays a change in its behavior at values of p for which xp = 1, as < f > 
becomes independent of p. 

In order to interpret these results, it shall be easier first to discuss on the kind of movement one 
would expect for cargo in the context of the considered ASEP. The mechanisms in ([3]) that define its 
elementary movements within each unit interval of time correspond to those of exchanging positions 
with a neighbor motor. The assigned hopping rates are such to promote, at a first moment, an 
accumulation of motors at one side (at the left side) of the cargo. Then cargo would be able to 
move backwards by exchanging position with these accumulated motors. By doing this, the motors 
end up transposed to the cargo's front. Because motors are assigned with an intrinsic dynamics - 
they move preferentially to the right - these motors at cargo's front will tend to disperse. Since 
the vesicle depends on such clusters of motors to develop a measurable velocity, it ends up moving 
mostly due to the motors accumulated at its back. 

The dependence of < f > on p in Fig. 2 confirms these expectations showing that < v > assumes 
only negative values, at all ranges of parameters. One could expect, in principle, that at high values 
of p there would be a balance between a tendency for maintenance of motors in front of the cargo, as 
excluded-volume become more important and eventually would be responsible for expressive motor 
accumulation at cargo's front. So, in principle, one could think that for sufficient high values of p 
there would be a chance for the cargo to develop a macroscopic movement towards the plus end 
of the microtubule, that is to the forward direction as well, < f > eventually displaying positive 
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values for p > x~^. According to the results, however, this does not happen either. The behavior 
of < V > predicted for high values of p can be understood by recognizing the formation of an 
"infinite" (macroscopic) cluster of motors behind the cargo with which it can always exchange 
positions. The formation of such infinite cluster would be a consequence of the phase transition 
(of condensation type) predicted for this system. The average velocity in this region of density 
becomes constant probably due to the uniformity of motor distribution along this cluster. 

We can then summarize these results by saying that the single cargo in this system of many 
motors develops, a backwards movement at any value of the hopping rates k,p and w, or density 
p. The magnitude of such velocity, however, is highly dependent on p and becomes constant at 
such values of p greater than a critical value p^ = x~^. It would be interesting then to test these 
predictions using data from experiment in vivo by monitoring the behavior of the motor density as 
cargo moves. 

Actually, data from Drosophila embryos [21], [2l] show that cargo velocity presents distinct 
behaviors depending on the stage of embryo development. It remains to investigate whether these 
changes could be associated to corresponding changes in the density of motors available at each 
of these stages. To our knowledge, there is limited information about the possible changes on 
the motor distribution along the microtubules due to the movement of cargo. The present study 
suggests that any investigation in this direction might be relevant to find ways to control cargo 
movement. 

Acknowledgments 

We acknowledge the financial support from Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo (FAPESP). 



14 



References 



[1] B. Derrida, M.R. Evans, The asymmetric exclusion model: exact results through a matrix ap- 
proach, Nonequilibrium statistical mechanics in one dimension. University Press, UK, Chapt. 
14, 277-304 (1997); B. Derrida, M.R. Evans, V. Hakim, V. Pasquier, An exact solution of a 
ID asymmetric exclusion model using a matrix formulation, J. Phys.A26, 1493-1517 (1993). 

[2] R.A. Blythe, M.R. Evans, Nonequilibrium steady states of matrix-product form: a solver's 
guide, J. Phys. A: Math Theor. 40, R333-R441 (2007). 

[3] D.Chowdhury, Traffic Flow of Interacting Self-driven particles: rails and trail, vehicles and 
vesicles. Physica Scripta T106, 13 - 18 (2003) . 

[4] A.Parmeggiani, T. Pranosch, and E. Prey, Phase Coexistence in driven one-dimensional trans- 
port. Phys. Rev. Lett. 90, 086601-1 - 086601-4 (2003). 

[5] A.Parmeggiani, T. Pranosch, and E. Prey, E. Totally Asymmetric Simple Exclusion Process 
with Langmuir Kinetics, Phys. Rev. E 70, 046101-1 - 046101-20 (2004). 

[6] S.Klumpp, R. Lipowsky, Traffic of Molecular Motors Through Tube-Like Compartments, J. 
Stat. Phys. 113, 233 -268 (2003). 

[7] O. Campas, Y. Kafri, K. B. Zeldovich, J. Casademunt, J.-F. Joanny, Collective Dynamics of 
Interacting Molecular Motors, Phys. Rev. Lett 97, 038101 (2006). 

[8] O. Campas, C. Leduc, P. Bassereau, J. Casademunt, J.-F. Joanny, J. Prost, Coordination of 
kinesin motors pulling on fluid membranes, Biophys. J. 94, 5009-5017 (2008). 

[9] J. Krug, P.A. Ferrari, J. Phys. A 29 (1996) L465-L471 (1996) . 

[10] M.R. Evans, Bose -Einstein condensation in disordered exclusion models and relation to traffic 
flow, Europhys. Lett. 36, 13-18 (1996). 

[11] A.Ploubidou, M. Way, Viral transport and the cytoskeleton, Curr. Opi. Cell Biol. 13, 97-105 
(2001). 



15 



[12] R.Kerstin, K. Dohner, B. Sodeik, Viral interactions with the cj^oskeleton: a hitchhiker's guide 
to the cell, Cell. Microb. 8, 387-400 (2006). 

[13] R.N. Cohen, M.J. Rashkin, X Wen, F.C. Jr. Szoka, Molecular motors as drug delivery vehicles. 
Drug Disc.Today: Technologies 2, 111-118 (2005). 

[14] D. G. Lichtenthaler, C. Goldman, Surfing at wave fronts: the bidirectional movement of cargo 
particles driven by molecular motors, preprint (2008) Cond Mat. arxiv.org/abs/0803.2732. 

[15] P. A. Ferrari, C. Kipnis, E. Saada, Microscopic Structure of TraveUing Waves in the Asymmetric 
Simple Exclusion Process, Ann. Prob. 19, 226 - 244 (1991). 

[16] T. Sasamoto, One-dimensional asymmetric simple exclusion process on a ring with a defect 
particle, Phys. Rev. E 61, 4980-4990 (2000). 

[17] H.-W. Lee, V. Popkov, D. Kim, Two-way traffic flow: exactly solvable model of traffic jam, J. 
Phys. A: Math. Gen.30, 8497-8513 (1997). 

[18] M. E. Fouladvand , H.-W. Lee, Exactly solvable two-way traffic model with ordered sequential 
update, Phys. Rev. E 60, 6464-6479 (1999). 

[19] B. Derrida, S.A. Janowsky, J.L. Lebowitz, E.R. Speer, Exact solution of the totally asymmetric 
simple exclusion process: shock profiles, J. Stat. Phys. 73, 813 - 842 (1993); 

[20] H. Kim, S.-C. Ling, G. C. Rogers, C. Kural, P.R. Selvin, S.L. Rogers, V.I. Gelfand, Microtubule 
binding by dynactin is required for microtubule organization but not cargo transport, J. Cell 
Biol. 176, 641-651 (2007). 

[21] M. A. Welte, S. P. Gross, M. Postner, S. M. Block, and E. F. Wieschaus, Developmental 
Regulation of Vesicle Transport in Drosophila Embryos: Forces and Kinetics, Cell 92, 547-557 
(1998). 

[22] D.H.U. Marchetti, P.A.F. da Veiga, T.R.Hurd, The 1/N-expansion as a perturbation about 
the mean field theory: a one-dimensional fermion model. Comm. Math. Phys. 179, 623-646 
(1996). 



[23] J. D. Murray, Asymptotic Analysis, Springer (1984). 

[24] S. P.Gross, Hilther and Yon: A Review of bidirectional Microtubule-Based Transport, Phys. 
BioLl, Rl-Rll (2004) . 



17 



Figure Caption 

• Fig.l - A configuration of the one-dimensional (discrete) ASEP model for interacting motors 
(gray) and cargo (black). Each particle, occupies a single site at each instant of time. The 
non-occupied (empty) sites are represented by line segments. The random processes are such 
that each motor is allowed to hop at rate k to its nearest neighbor on the right if it is empty. 
The motors and cargo can exchange places at a rate p (motor jumping to nearest neighbor 
at the right ) or q (motor jumping to the nearest neighbor at the left). 

• Fig. 2 - Average velocity of cargo as a function of motor density p, at various values of param- 
eter p, for fixed w — 3 and A; = 1, as indicated. The predicted phase transition is illustrated 
by the change in the behavior oi < v > that assume a constant value ior p > — 1/x. 
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